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Introduction 

This report covers the author's research activities at the NASA Langley 
Research Center during the ten-week period June 11 to August 17, 1973. This 
work was supported under NASA Contract No. NASl-1 1 707-21 , and was carried out 
in the Hypervelocity Impulse Facilities Section of the Advanced Entry Analysis 
Branch, Space Systems Division. Mr. Jim Jones served as Immediate supervisor. 
The effort involved work In the three areas described briefly below. 

Correlation of Aerodynamic Heating Rates 

This work Involved a theoretical study of heat transfer for zero pressure 
gradient hypersonic laminar boundary layers for various gases with particular 
application to the flows produced In Langley expansion tube facility, A 
simple correlation based on results obtained from solutions to the governing 
equations for five gases was formulated. See the attached first draft of a 
paper entitled, "Correlation of Heat Transfer for the Zero Pressure Gradient 
Hypersonic Laminar Boundary Layer for Several Gases," which will be submitted 
to the AIAA Journal. 

Heat Transfer Measurement In Short-Duration Flows 

The author served as a consultant on data reduction techniques in relation 
to the measurement of heat transfer rates in the expansion tube facility, 
techniques in use by the Langley scientists which were previously developed 
^y the author for reduction of data obtained from thin-fllm resistance thermo- 
meter heat-flux gages were extended for use In conjunction with th1n-film 
thermocouple gages. 

^ck Tube Boundary Layer Analysis 

The author served as a consultant on the subject of high temperature 
^3fn1nar boundary layer flows in shock tubes. Particular attention was directed 
toward the laminar boundary layer on shock tube splitter plates In carbon 
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Hoxide flows generated by high speed shock waves. Computer analysis of the 
;plitter plate boundary layer flow provided information that is useful in 
nterpreting experimental data obtained in shock tube gas radiation studies. 
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CORRELATION OF HEAT TRANSFER FOR THE ZERO PRESSURE GRADIENT 
HYPERSONIC LAMINAR BOUNDARY LAYER FOR SEVERAL GASES 

By William J. Cook* 

NWENCLATURE 

Chapman-Rubesin coefficinet. pp/p p 

6 6 

Static enthalpy, J/kg 

Recovery enthalpy, h + ruV2 

e e 

Dimensionless enthalpy, (h - h )/h 

e e 

Thermal conductivity, W/(m K) 

Pressure, atm 

Heat transfer rate per unit area, (heat flux), W/m2 
Recovery factor 

Reynolds number, xu o /u 

e^e'^e 

Mach number 
Temperature, K 

Exponent: cone, s = 1; wedge, plate, s = 0 
Stanton number, -q^PgUgCh^ - h^) 

X component of velocity, m/s 
y component of velocity, m/s 
Space coordinate parallel to wall, m 

Space coordinate perpendicular to wall, m 
u/u. 
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Y Specific heat ratio, Cp/c^ 

X Dimensionless shear stress, y(8u/ay)/p^u| 

Viscosity coefficient, kg/ (m s) 
p Density, kg/m^ 

^ Shear stress variable, Eq. (3) 

a Effective Prandtl number 

Subscripts 

e Boundary layer edge 

w Wall 

INTRODUCTION 

Current research in the area of convective heating in high speed flows 
includes studies using gases and gas mixtures to simulate aerodynamic heating 
in planetary atmospheres. Test facilities such as the shock tunnel and the 
expansion tube are useful for generating the required high-*energy flows for 
investigating aerodynamic heating over a range of high Mach numbers. Simple 
configurations such as the sharp flat plate, the wedge, and the sharp cone 
that avoid boundary layer pressure gradients are frequently used as basic 
body shapes in such studies. This note is concerned with correlation of 
hypersonic aerodynamic heating parameters for these body shapes for several 
gases under test conditions encountered in these facilities and for conditions 
under which the boundary layer is laminar. 

Testing in the above-mentioned facilities characteristically involves 
models with cold surfaces (relative to the recovery temperature) since, due 
to short testing times involved, little change from the initial room temperature 
of the model surface occurs. Hence, considerations here will be directed mainly 
toward flows with the cold wall condition. Further, only the boundary layer 
flow downstream of the region of Interaction of any leading edge shock wave 
2nd the boundary layer that occurs near the leading edge will be treated. 

To a good approximation the self-similar laminar boundary layer equations 
Tor flows with zero pressure gradient are then applicable. 
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In order to assess the influence of the characteristics of various real 
gases, hypersonic heating rates in air, argon, carbon dioxide, helium, and 
hydrogen have been theoretically studied. These gases as ideal gases range 
in molecular structure from monatomic to tri atomic, the molecular weights 
range from 2 to 44 and values of y vary from 1.28 to 1.67. In high Mach 
number flows real gas effects can play an important role since, because of 
pronounced dissipation effects, the temperature in the boundary layer can 
reach large values. Real gas effects in hypersonic flows may be present 
to varying degrees depending on the gas considered, e.g., carbon dioxide 
dissociates at relatively low temperatures while helium exhibits ideal 
gas behavior over a wide range of temperature. In view of these facts and 
the range of gas characteristics considered, there is no assurance that a 
precise universal correlation for heat transfer exists for a wide range of 
flow conditions. However, insight for the present problem can be gained 
by examining the boundary layer equations and solutions thereto for various 
gases. 


ANALYSIS 


The method employed here to solve the governing boundary layer equations 
is that of Crocco (e.g., see Ref. 1), in which the boundary layer equations 
in terms of x and y as the independent variables are transformed to 
equations that have x and u as the independent variables. For real gas 
flows in thermochemical equilibrium, the laminar boundary layer equations 
with 3 = u/u and x as the independent variables and x = v(3u/9y)/p u^ 

6 G € 

and H = (h - h.)/h^ as dependent variables are: 

c c 

Momentum, 


+ B 


“ Py__l 

t 

pyB 

e e . 

+ 



x 



= 0 


Energy, 




L<J Bi p|v 


( 1 ) 


( 2 ) 
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where s = 1 for the case of the cone and s “ 0 for the plate and the wedge. 
In Eq. (2) a = a{h,p) is the effective Prandtl number that serves to incorp- 
orate real gas effects with regard to energy transport. Continuity is 
implicitly satisfied in Eqs. (1) and (2). The boundary conditions for these 
equations are, 

3t/30 = 0, H = H . at 0 = 0 

W 

T = 0, H f 0, at 0 = 1 

The boundary condition on t at 0 = 0 results from u = v = 0 at the wall. 

By employing a separation of variables solution of the form t(x,0) = 0(b)/X(x), 
the shear stress becomes, 

t(x, 0) = 0(0)/{Re)^^^ (3) 

After Crocco, the term in Eq. (2) is taken as zero. Substituting Eq. (3) 
into Eqs. (1) and (2), 

+ 3®Cg/2 =0 (4) 

%% - (5) 

where C = py/p y . The boundary conditions become 

= 0, H = at 0 = 0 
0 = 0, H = 0, at 0 = 1 

In integrated form Eqs. (4) and (5) can be written as. 



( 6 ) 
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in U(X)] 


°x ‘’x 


where t(g) = [0(6)]^”'’^^^. Equations (6) and (7) are coupled since C = C(h,p). 

The heat flux q = -k(3T/ay) for the gas becomes 
w w 


’w 


= “0, h^/a,, (Re) 

w ' g'w e e e' w ' ^ 


1/2 


( 8 ) 


where 0^ and (H^)^ are obtained from simultaneous solution of Eqs. (6) and 
(7). Substitution of Eq. (8) Into the Stanton number expression yields. 


St(Re) 


1/2 . 






1/2 u2 


2h 


- «w' 


(9) 




1/2 


where the recovery factor r has been taken as 

It is useful to examine the equations that lead to the terms in Eq. (9) 
in order to determine how the various flow quantities and gas properties 


influence St (Re) 


1/2 


To solve Eqs, (6) and (7) for 0^ and 




certain 


quantities must be specified. For a given gas the variations of the properties 
a, T, and py with h at the specified pressure are required. Specification 


of T and 
w 

c = PP/PgUe 


T then establish h and h and therefore H . 
e we w 

and cr can then be determined as functions of H. 


In addition. 
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Finally, with specification of u^, the term u^/h^ in Eq. (7) is 
f. e e e 

[established. Thus, if T , T , u , and the gas properties as a function 
of h are specified, sufficient information is available to solve Eqs. (6) 
and (7) for 0^ and Therefore, if T^, T^, u^, and the pressure 

p are specified, and various gases are considered, any variation that occurs 
in 0,^ and (H ) will arise from differences in the relations for T, py, 
and 0 vs h for the various gases. It should be noted however that fixed 
values for T and T do not yield the equal values of for all gases 
except for the case T^ = T . Therefore, if comparisons are to be made between 
the solutions to Eqs. (6) and (7) for various gases at conditions other than 
T^ = Tg, they should be made at the same values of T^ and rather than 
at the same T^ and T^. On this basis, specification of fixed values of 

P» would permit the variations of py and a vs h 

on 0^ and (Hq)., for various gases to be assessed. The terms u|/h_ and 

W p W 6 6 

also appear in the denominator of Eq. (9). The term u£/h^ is related 

W 6 6 

to the Mach number. For ideal gases u£/h = M2(y - 1). Thus for large 

Mach numbers i^^e Prandtl number in the denominator of 

Eq. (9) for the cold walls considered here is related to y« (According to 

2 

the theory of Eucken the Prandtl number for gases can be closely approximated 
by 4y/(9y - 5).) Therefore, for high Mach numbers the denominator of Eq. (9) 
is related to the Mach number and y- Computations performed to examine the 
influence of variations of gas properties on Eq. (9) and the terms therein 
are presented in the next section. 

RESULTS 

Figure la shows solutions to Eqs. (6) and (7) for a plate (s = o) in terms 
of t^(g) and H(0) obtained for carbon dioxide and helium for the values of 
^w* P* ^w figure. The numerical method employed 

to solve Eqs. (6) and (7) is outlined in ref. 3. Figure lb shows C(g) and a($) 
for each of the gases. Property values for carbon dioxide were obtained from 
refs. 4 and 5. Thermodynamic properties for helium were computed from ideal 
9as equations. Also the relation y and a constant Prandtl number 

(o = 0.688) were employed for helium. 

It is noted that the influence of the difference between the curves for 
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C{e) and <j(a) on the profiles for <^( 3 ) and H(e) is small, but that 

the values for noted in the figure differ by about 20 percent, 

W p w 1 

However, the values of St(Re) ' computed by Eq. (9) for the two gases 
differ by less than one percent. This suggests that St (Re) ' and u|/h^ 
are the correlating variables. Figure 2 presents results in terms of St(Rg 
vs u£/h. for air, argon, carbon dioxide, helium, and hydrogen. The term 
permits the results to be applied to the cone or the plate and wedge. 

The properties for argon and hydrogen were obtained from ref. 4. Those for 
air were taken from ref. 6, The Mach number and u^/h^ ranges applicable for 
each gas are shown in the figure. Results obtained for H., = 0 and H = -0.68 
are indicated by the open and filled symbols respectively. These values of 
correspond to T^/T^ = 0 and T^/T^ * 3, and cover the range of T^/T^ 
commonly encountered in shock tunnel and expansion tube flows. With the 
exception of carbon dioxide the results are relatively insensitive to the 
value of H . The results in Fig. 2 were obtained for p = 0.01 atm. However, 

W 

computations at p = 1 atm for carbon dioxide (for which variation in pressure 
was expected to have the most influence) yielded results that are within one 
percent of those shown in Fig. 2, indicating that there is no significant 
influence of pressure for the conditions considered here. 

For simplicity a straight line has been fitted through the points in 
Fig. 2. The computed results deviate by at most five percent from this 
curve. Thus for most engineering purposes this simple relation should be 
sufficiently accurate for the gases considered for Fig. 2. In view of the 
wide range of gas characteristics considered here it is reasonable to 
expect this relation would be applicable with the same accuracy to other 
gases and gas mixtures provided the cold wall condition is maintained and 
is within the range noted in Fig. 2. 
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